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We present a detailed analysis of the fluctuation dissipation theorem (FDT) close to the glass transition in 
colloidal suspensions under steady shear using mode coupling approximations. Starting point is the many- 
particle Smoluchowski equation. Under shear, detailed balance is broken and the response functions in the 
stationary state are smaller at long times than estimated from the equilibrium FDT. An asymptotically constant 
relation connects response and fluctuations during the shear driven decay, restoring the form of the FDT with, 
however, a ratio different from the equilibrium one. At short times, the equilibrium FDT holds. We follow 
two independent approaches whose results are in qualitative agreement. To discuss the derived fluctuation 
dissipation ratios, we show an exact reformulation of the susceptibility which contains not the full Smoluchowski 
operator as in equilibrium, but only its well defined Hermitian part. This Hermitian part can be interpreted as 
governing the dynamics in the frame comoving with the probability current. We present a simple toy model 
which illustrates the FDT violation in the sheared colloidal system. 

PACS numbers: 82.70.Dd, 64.70.P-, 05.70.Ln, 83.60.Df 



I. INTRODUCTION 

The thermal fluctuations of a system in equilibrium are di- 
rectly connected to the system's response to a small external 
force. This connection, manifested in the fluctuation dissipa- 
tion theorem (FDT), lies at the heart of linear response theory. 
The FDT in equilibrium connects the correlator (t) with 

the response function, the susceptibility Xfg {t) (both defined 
below), and reads 
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Eq. ((TJ states that the relaxation of a small fluctuation is in- 
dependent of the origin of this fluctuation: Whether induced 
by a small external force or developed spontaneously by ther- 
mal fluctuations, the relaxation of the fluctuation cannot dis- 
tinguish these cases. 

The most famous example for the FDT is the Einstein re- 
lation connecting the diffusivity of a Brownian particle to its 
mobility |Q]]. The FDT is of importance for various appli- 
cations in the field of material sciences since for example 
transport coefficients can be related to equilibrium quantities, 
i.e., the fluctuations of the corresponding variables 0. It 
was first formulated by Nyquist in 1928 yfl as the connec- 
tion between thermal fluctuations of the charges in a conduc- 
tor (mean square voltage) and the conductivity. 

In non-equilibrium systems, this connection is not valid in 
general and much work is devoted to understanding the rela- 
tion between fluctuation and response functions. This rela- 
tion is often characterized by the fluctuation dissipation ratio 
(FDR) Xfg(t) defined as 



Xfg(*) 



k R T dt 



C fg (t) . 
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Close to equilibrium, one recovers the FDT in Eq. ([l} with 
1. In non-equilibrium, Xf g (t) deviates from unity. 



This is related to the existence of non-vanishing probabil- 
ity currents (see Eq. (fl~5T > below); FDRs are hence consid- 



ered a possibility to quantify the currents and to signal non- 
equilibrium [0]. The violation of the equilibrium FDT has 
been studied for different systems before as we want to sum- 
marize briefly. 

The general linear response susceptibility for non- 
equilibrium states with Fokker-Planck dynamics [j]!] has been 
derived by Agarwal in 1972 JfJ. It will serve as exact starting 
point of our analysis (see Eq. (foi l below). The susceptibil- 
ity is given in terms of microscopic quantities, which cannot 
easily be identified with a measurable function in general in 
contrast to the equilibrium case. For a single driven Brownian 
particle (colloid) in a periodic potential, the FDT violation for 
the velocity correlation has been studied in Ref. 0. There it 
was possible to compare the microscopic expressions success- 
fully to the experimental realization of the system 0]. 

Colloidal dispersions at high densities exhibit slow cooper- 
ative dynamics and form glasses. These metastable soft solids 
can be easily driven into stationary states far from equilibrium 
by already modest flow rates. Spin-glasses driven by non- 
conservative forces were predicted to exhibit nontrivial FDRs 
in mean field models @]. It is found that at long times the 
equilibrium form of the FDT (Eq. (Q]i) holds with the temper- 
ature T replaced by a different value denoted effective tem- 
perature T e ff, 



-1 8 

k B T cS dt 



C{i). 



(3) 



This corresponds to a time independent FDR Xf(t) = Xf 
at long times during the final decay process, where t is time 
rescaled by the timescale of the external driving. In detailed 
computer simulations of binary Lennard- Jones mixtures by 
Berthier and Barrat 11211 . this restoration of the equi- 

librium FDT was indeed observed: For long times, the FDR 
Xf is independent of time. Its value was also very similar 
for the different investigated observables, i.e., Xf = X = 
T/T e s is proposed to be a universal number describing the 
non-equilibrium state. T e g was found to be larger than the 
real temperature, which translates into an FDR smaller than 
unity. Further simulations with shear also saw T e g > T 
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lfl3l [Til IT5I [l6ll . but the variable dependence was not stud- 
ied in as much detail as in Ref. fioll . and partially other def- 
initions of T c ff were used. In Refs. H Q3 [nil it is argued 
that T c ff agrees with the effective temperature connected with 
the FDT violation in the corresponding aging system lfl8i[l9tl . 
This has not yet been demonstrated for different temperatures. 
Note that the system under shear is always ergodic and aging 
effects are absent. The fluctuation dissipation relation of aging 
systems using mode coupling techniques was investigated in 
Ref. 12011 . Recently T e g was also connected to barrier crossing 
rates Il2ll]r eplacing the real temperature in Kramers' escape 
problem [50. A theoretical approach for the effective tempera- 
ture under shear in the so called "shear-transformation-zone" 
(STZ) model is proposed in Ref. lf22ll . Different techniques 
(with different findings) to measure FDRs in aging colloidal 
glasses were used in Refs. 1I23I I24II25I1 . No experimental re- 
alization of an FDT study of colloidal dispersions under shear 
is known to us. An overview over the research situation (in 
2003) can be found in Ref. @]. 

Interesting universal FDRs were found in different spin 
models under coarsening lEil I27L EM I29I l30ll and under shear 
ll3lll . At the critical temperature, a universal value of X = | 
has been found e.g. in the n-vector-model for spatial dimen- 
sion d > 4. In d = 3, corrections to this value can depen- 
dent on the considered observable 13211 . See Ref. lO for an 
overview. Yet, the situation for structural glasses has not been 
clarified. Also, the connection between structural glasses, spin 
glasses and critical systems is unclear. 

In this paper, we present the study of the violation of 
the equilibrium FDT for dense colloidal suspensions under 
shear. It is a comprehensive extension of our recent Letter 
on the same topic J3, but also provides a number of new 
results and discussions. We build on the MCT-ITT approach 
JUl H2I HH] (reviewed recently J39D) based on mode cou- 
pling theory. This approach allows us to derive quantities 
which are directly measured in experiments and simulations 
ll40l RH l42ll and the properties of specific observables can be 
described. We will hence be able to study the non-equilibrium 
FDT for different observables as measured in simulations, and 
possible differences for different variables can be detected. In 
the main text, we will follow the calculation as presented in 
Ref. l34ll in detail. It leads to a time independent FDR Xf dur- 
ing the whole final relaxation process whose value is universal 
in the simplest approximation, Xf = |. We will also derive 
corrections to this value which depend on the considered ob- 
servable. In Appendix|A] we will additionally show a different 
analysis of the extra term in the FDT following more standard 
routes of MCT and projection operator formalisms. It is in 
qualitative agreement with the results shown in the main text 
and it allows us to estimate the size of the correction terms 
which are neglected in the main text and to see that they are 
small. 

The paper is organized as follows. In Sec. [Til we will in- 
troduce the microscopic starting point and give the definitions 
of the different time dependent correlation and response func- 
tions. In Sec. Hill we will introduce the different contribu- 
tions to the non-equilibrium term A\f g (t) in the susceptibil- 
ity. These different contributions are approximated in Sec.lIVI 



The approximations for the time dependent correlation func- 
tions will be shown in Sec. [V] In Sec. |VI| we will present 
the final extended FDT connecting the susceptibility to mea- 
surable correlation functions and discuss the FDR as function 
of different parameters. In Sec. I VIII we show an exact form 
of the susceptibility which involves the Hermitian part of the 
Smoluchowski operator and the restoration of the equilibrium 
FDT in the frame comoving with the probability current. The 
final discussion, supported by the FDR analysis in a simple toy 
model will finally be presented in Sec. IVIIII In Appendix lAl 
we derive the expressions for the susceptibility in an approach 
based on the Zwanzig-Mori projection operator formalism. 



n. MICROSCOPIC STARTING POINT 

We consider a system of N spherical Brownian particles of 
diameter d, dispersed in a solvent. The system has volume 
V. The particles have bare diffusion constants Do- The inter- 
particle force acting on particle i (i = 1 ... N) at position r ; 
is given by Fj = — d/driU({rj}), where U is the total po- 
tential energy. We neglect hydrodynamic interactions to keep 
the description as simple as possible. These are also absent in 
the computer simulations iflOfl to which we will compare our 
results. 

The external driving, viz. the shear, acts on the particles 
via the solvent flow velocity v(r) = 72/x, i.e., the flow points 
in x-direction and varies in y-direction. 7 is the shear rate. 
The particle distribution function ^(T = {rj},i) obeys the 
Smoluchowski equation l37ll43ll . 

$*(r,t) = n*(r,t), 

si = n e + sn = % ■ [9i - Fi - « ■ T t ] , (4) 

i 

with k = 7xy for the case of simple shear. SI is called the 
Smoluchowski operator (SO) and it is built up by the equi- 
librium SO, Sl e = J2i 9i ■ [di — Fj] of the system with- 
out shear and the shear term Sil — — J2i 9% ■ K • r,. We 
introduced dimensionless units for space, energy and time, 
d = fcgT = Do = 1. The formal H-theorem ^ states 
that the system reaches the equilibrium distribution *f? e , i.e., 
Oe^e = 0, without shear. Under shear, the system reaches 
the stationary distribution & s with Qty s = 0. Ensemble av- 
erages in equilibrium and in the stationary state are denoted 

(...)= J dT¥ e (r)..., (5a) 
= J dT*.(r)..., (5b) 

respectively. In the stationary state, the distribution function 
is constant but the system is not in thermal equilibrium due to 
the non-vanishing probability current l37tl . 

j? = [-di + F» + k ■ ri ]* s = j&s. (6) 
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FIG. 1: Definition of the waiting time t w and the correlation time t 
after switch-on of the rheometer. 



A. Correlation functions 



B. Response Functions 

The susceptibility \f g describes the linear response of the 
stationary system to an external perturbation. Note that the 
term 'linear response' does not correspond to the shear, but 
to the additional small test force h e (t) acting on the particles. 
Because the system is always ergodic due to shearing, the lin- 
ear response will always exist in contrast to un-sheared glasses 
l44l,l45ll . where a finite force is needed to mobilize the parti- 
cles. Formally, the susceptibility Xf g (i) describes the linear 
response of the stationary expectation value of g to the exter- 
nal perturbation h e (t) which shifts the internal energy U to 

U-r he(t), 



Dynamical properties of the system are probed by time de- 
pendent correlation functions. The correlation of the fluctua- 
tion Sf = /— of a function /({r$}) with the fluctuation 
of a function g({ri}) is in the stationary state given by l37ll 



(7) 



(7) 



Here, = + F * + r « ' kT ] ' ®i is the adjoint SO that 

arose from partial integrations. Cf s (t) is called the station- 
ary correlator, it is the correlation function which is mostly 
considered in experiments and simulations of sheared suspen- 
sions. At this point, we would like to introduce three more 
correlation functions which will appear in this paper. The 
transient correlator ci „ is observed when the external shear 



(8) 



is switched on at t = IB 



It probes the dynamics in the transition from equilibrium to 
steady state and is the central object of the MCT-ITT approach 
11361. l38ll . In the general case, where the correlation is started 
a period t w , namely the waiting time, after the rheometer was 
switched on, one observes the two-time correlator Cf g (t, t w ), 
see Fig.Q] 

Q s (t,t w ) = Cg 5 (t)+7 f " ds (<J xy e n1s Sre nU Sg) . (9) 



Eq. (0 follows with the waiting time dependent distribution 
function JH, 



dse ns m c 



(10) 



with partial integrations when averaging with ^(T,t w ). For 
t w = 0, the two-time correlator equals the transient correlator. 
For long waiting times, it reaches the stationary correlator, 
Cf s (t,t w — > 00) — » Cfg(i), and Eq. (0 becomes the ITT 
expression for C{ s (t) ll37ll . Without shear finally, one observes 
the equilibrium correlation, 



a 



(e) 

fg 



Sf 



(ID 



(g) w le) (t) {g) w = / dt'xfeC* - t')h e {t>) + o(hl) 



(12) 

To derive the microscopic form of the susceptibility, one con- 
siders the change of the stationary distribution function s 
under the external perturbation. One finds f5l, llil, l37tl 




(7) 



(13) 



If one replaces *f> s by ^ e and Qt by fit m £q. ([T3l , the equi- 
librium FDT in Eq. (Q]i follows by partial integrations, 



f*Qte^g 



(14) 



In the considered non-equilibrium system, where detailed bal- 
ance is broken and the nonzero stationary probability current 
in Eq. © exists, the equilibrium FDT ([T]i is extended as we 
see now. The above expression, Eq. (1131 , can be rewritten 
(with adjoint current operator j\ = di + F, + r-j • k t ) to 



Axfg(i) = Xf g (t) + C-fe(t) = - ( J2 St ■ °-£- ^ 



(7) 



(15) 

Note that the new term in the FDT, A%f g (t), in the following 
called violating term, is directly proportional to the station- 
ary probability current. A deviation of the fluctuation dissipa- 
tion ratio, Eq. (O, from unity, the value close to equilibrium, 
arises. 

The extended FDT in Eq. ( fT~5l > has been known since the 
work of Agarwal We will analyze it for driven metastable 
(glassy) states and show that the additive correction Axf g (£) 
121 |H |46[] leads to the nontrivial constant FDR at long times, 
as was found in the simulations. One can always express the 
FDT violation in terms of an additive as well as a multiplica- 
tive correction. The nontrivial statement for driven metastable 
glasses is that the multiplicative correction is possible with a 
time independent factor at long times. Specifically, we will 
look at autocorrelations, g = f of functions without explicit 
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advection, / = f({yi, z,}), where the flow-term in the current 
operator j\ in < fT3T > vanishes. For variables depending on Xi, 
the equilibrium FDT is already violated for low colloid densi- 
ties as seen from the Einstein relation: The mean squared dis- 
placement grows cubically in time J47ll (Taylor dispersion), 
while the mobility of the particle is constant. 

In contrast to the equilibrium distribution, oc e~ u , the 
stationary distribution is not known and stationary averages 
are calculated in the integration through transients approach 
(ITT) JH (compare Eq. ([10>), 



5f*5tfe uH 8f) = -^C t (t,t w ) 



(19) 



(...)« = (■• ->+7 



) 



ITT simplifies the following analysis because averages can 
now be evaluated in equilibrium, while otherwise non- 
equilibrium forces would be required l48h . E.g. due to 
di^ e = Fi^e, the expression (fT~5T > vanishes in the equilib- 
rium average and reduces to 



where from now on, we consider fluctuations from equilib- 
rium, Sf = f — {/). The constant (/) cancels in (fT~8b . The 
second equal sign in Eq. $1% follows with partial integrations 
(recall Sfl^e = a xy ^ e and [<5fit,<5/] = 0). The intriguing 
connection to the waiting time derivative follows by compar- 
ison of the right hand side of the first line of Eq. (fT9l > with 
Eq. ©. 

The second term in Eq. ( fTST l describes the time derivative 
of the difference between stationary and transient correlator, 
compare Eq. (0 with t w — > oo. The last term in Eq. ( TT~8b 
has no physical interpretation and we denote it by A\ { 
hence have 



(3) 



We 



t w =0 



Ci{t)-C\ l> {t) + 



/3) 



(20) 



Axf(*) 



ds 



fs J2( 9 * + F ') • JT^** ) ■ where A X f ) = / °° ds(a xy e n1 s (fit f*)e nU /). In the fol- 



>xy* 



dr 



(16) 

Eq. ( TToT ) is still exact and we will in the following develop ap- 
proximations for it. In the main text, we will follow the deriva- 
tion as presented in Ref. l34ll . In Appendix lAl we present an 
alternative derivation with Zwanzig Mori projections. We will 
show in Sec. IA4I that the two approaches are in qualitative 
agreement. 



III. THE VIOLATING TERM 

We want to analyze the violating term A%f (i) in more de- 
tail. It can be split up into terms containing the Smoluchowski 
operator instead of the unfamiliar operator jt. This can be 
done with the following identity for general functions /({r^}) 
andg({r l }), 

£j1 ■ §^.9 = \ [nVs - r^g + (fitr ). 9 ] . a?) 



If we apply this identity to Eq. ( fT6] l with g = e nt */, we get 
the following three terms, 



lowing subsections, we will look at the different terms more 
closely. 



IV. APPROXIMATIONS FOR THE VIOLATING TERM 

A. The waiting time derivative 

In order to approximate the waiting time derivative in 
Eq. ( fT9b , we note its connection to time derivatives of cor- 
relation functions, 



Sf* c5fi f e n t 6f 



C?\t)-Urn\e nU 5f 



(21) 



Ax f (i) 



/ ds{a xy e n ' s pf* - + (fit/*)]e° Tt /). 



(18) 

The first term in Eq. ( TT~8b contains a derivative with respect 
to s and the s-integration can immediately be done. We find 
that the first term in Eq. dT8l (without the factor h) exactly 
describes the derivative of Cf (t,t w ) with respect to t w at t w = 
0, 

- 7 J ds (a xy e n ' s tf5f*e nU 8f) = 7 (<J xy Sre nh 6f 



The time derivative of the transient correlator cjp (t) is split 
into two terms, one containing the equilibrium operator fit, 
the other one containing the shear term <5fit. We will reason 
the following: The term containing fit is the derivative of the 
short time, shear independent dynamics of the transient cor- 
relator down on the plateau (compare Fig. [3]below), i.e., the 
derivative of the dynamics governed by the equilibrium SO 
fi e . The term containing 5W, i.e., the waiting time derivative, 
follows then as the time derivative with respect to the shear 
governed decay from the plateau down to zero. 

The equilibrium derivative fi£<5/* in the last term of d2Tb 
de-correlates quickly as the particles loose memory of their 
initial motion even without shear. In this case, the latter term 
is the time derivative of the equilibrium correlator, Cf C \t). 
A shear flow switched on at t = should make the particles 
forget their initial motion even faster, prompting us to use the 



ntt c o+t 



approximation e nTt w e n ='Pfe a "="e 
SDiSpSf)- 1 ^/. We then find 



, with projector P{ 



(Sf*e- n - t e nh Sf) 



(22) 
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The first average on the right hand side is the time derivative 
of Cf (i). The second average is not known. Applying the 
same approximation e fi+t 
correlator, we have 



s nl*p fe - 



n e t e n+i t0 ^ transient 



~(<r#'ti y (5f , Sf) 1 



(23) 



Combining the two equations, we find for the last term in 
Eq. (ED 



cf ft) 



(24) 



This term is then assured to decay faster than without shear. 
Now we can give the final formula for the waiting time deriva- 
tive, 



d 

-z—C{(t,t w ) 
ot w 



u\ C i ft) 



Cf>(t)-Cp(t) 



(25) 



This is our central approximation whose consequences for the 
FDR will be worked out in Sec. [VI] The quality of approxima- 
tion d25l l has recently been studied in detailed simulations, and 
qualitative and quantitative agreement was found for two dif- 
ferent simulated super-cooled liquids Ir49ll . As argued above, 
the last term in $25[ will be identified as short time derivative 
of Cf , connected with the shear independent decay, where 
the transient correlator equals the equilibrium correlator. Con- 
sequently, gf-Cf (t, t w )\t m =o will turn out to be the long time 

derivative of Cf, connected with the final shear driven de- 
cay. This captures the additional dissipation provided by the 
coupling to the stationary probability current in Eq. (fTBT l. The 
approximation in Eq. $25[ is also reasonable comparing it to 
the expected properties of the waiting time derivative: For 
long times, t — > oo and 7^0 with jt = 0(1), one has 

0f ^ ft) — in glassy states and the waiting time derivative is 
equal to the time derivative of the transient coiTelator. Vary- 
ing the waiting time or the correlation time has then the same 
effect on the transient correlator. It is for small waiting times 
a function of j(t + t w ) since (t + t w ) measures the time since 
switch-on 113011 . 



B. The other terms in Eq. QUl 

The second term in Eq. ( TT~8T > has a physical interpretation 
as well: It is the time derivative of the difference between 
stationary and transient correlator, see Eq. (O, 
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ds(a 



.ry< 



s 5f*rte nJt Sf) = C f (t) - C [ p(t). (26) 



(3) 

The last term, A\\ , has yet no physical interpretation. At 
t = 0, it cancels with the second term. It is a demanding task 
to estimate the contribution of the different terms to Axi(t). 
This can be done in an MCT analysis for density fluctuations 





0.2 r 




0.1 








II 


-0.1 






D" 


-0.2 






-0.3 - 




-0.4 - 




-0.5 



• 




□ □ □ □ 




r .snLaf />' 

/ H m * a /' 
3" / ' 



10 

q 



15 



20 



FIG. 2: The estimates for the Laplace transform Ax q (z = 0) of 
the violating term in density susceptibilities for q = qy. We show 
the coherent (incoherent) case with data points connected (uncon- 
nected) by lines. Solid squares show the estimate for the first term, 
the waiting time derivative, solid spheres show the sum of the two 
other terms. Open squares show the sum of all three. The solid 
(dashed) line shows the prediction of Eq. j25b for the waiting time 
derivative for the coherent (incoherent) case; compare with the solid 
squares and see main text. The picture is very similar for the z- 
direction 15011. 



as presented in AppendixlAl We want to briefly summarize the 
results for the contributions of the different terms as found in 
AppendixlAl For coherent, i.e., collective density fluctuations, 
we have / = £> q = £\ e tq ri pill . For incoherent, ie., single 
particle fluctuations, one has / = f?q = e' q r = with r s the po- 
sition of the tagged particle. We denote all normalized density 
functions with subscript q, the normalized transient density 
correlator is denoted $ n 



3611 . We find that the violating 
For long times, jt = 0(1), we 

as 

is shown in Appendix|A] We find 



term is zero for jt <C 1 

can estimate the different contributions in terms of < & q (t) 



l -^-C(tt 
2 dt w Ccl[t ' U 



t w =o 



fiW$ q (t), (27a) 



C q (i) - * q (i)l +A Xq 3) W » n q 2+3) $ q (t),(27b) 



where the functions n q oc I7I for small shear rates in 
glassy states. In Fig. [2] we show the Laplace transform 
£,{Axq(t)}(z) = J °° dte- zt A X ^it) at z = for the differ- 
ent contributions to Ax q as estimated in AppendixlAl We see 
that, according to the estimates in AppendixlAl the first term 
in Eq. <fl~8T > is the dominant contribution to the violating term. 
It is larger than the two other terms, which additionally par- 
tially cancel each other. We will in the main text neglect the 
sum of second and third term, this will give good agreement 
to the data in Ref. fioll . In Fig. [2] we also show the prediction 
of Eq. d25l l for the waiting time derivative, which is at z = 
simply given by minus the height of the glassy plateau, see 
Eq. (|A17| >. We see that our estimate in Appendix [A] agrees 
qualitatively and also semi-quantitatively with Eq. 
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It appears reasonable to conclude that the waiting time 
derivative is larger than e.g. the second term ( |26] l, since it is 
equal to the time derivative of the transient correlator at long 
times, whereas (l26l i is the difference of two very similar func- 
tions. It is equal to the third term at t = which lets us expect 
that also Ax^ 3 ' is small. 



Eq. (f30b is a short version which neglects the waiting time 
dependence of the t = value of the unnormalized two-time 
correlator. An extended version including this effect can be 
formulated j49[ l50ll but it is more involved and would change 
our results only marginally. Thus, we continue with Eq. (130b 
where the first factor on the right hand side is the normalized 
integrated shear modulus 



V. APPROXIMATIONS FOR CORRELATION FUNCTIONS 

A. ITT equations for transient correlators 

The known ITT solutions for the transient correlators will 
be the central input for our FDR analysis. In order to visualize 
its time dependence, we will use the schematic -model 
of ITT, which has repeatedly been used to investigate the dy- 
namics of quiescent and sheared dispersions [36], and which 
provides excellent fits to the flow curves from large scale sim- 
ulations ll53ll . It provides a normalized transient correlator 
C M = ^K*)> as we H as a quiescent one, representing co- 
herent, i.e., collective density fluctuations. The equation of 
motion reads 13611. 



= $(i) + r + J dt'm('y,t j , (28a) 

t} = T+W ^ + 2Ale)m + v ^ 2{t) } ' (28b) 

with initial decay rate T. We use the much studied values 
i>2 = 2, vf = v^iy^Mi ~ 1) with the glass transition at 
e = 0, and take m(0, t) in order to calculate quiescent (7 = 0) 
correlators llsill . Positive values of the separation parameter e 
correspond to glassy states, negative values to liquid states. 

In order to study the g-dependence of our results, we will 
use the isotropic approximation l36ll for the normalized tran- 
sient density correlator. For glassy states, the final decay from 
the glassy plateau of height f q is approximated as exponential, 



(29) 



where the amplitude h q is also derived within quiescent MCT 

m. 



B. Two-time and stationary correlator 

We will need to know the difference between stationary and 
transient correlators in order to be able to study the FDR in 
detail. Here we derive an approximate expression for the two- 
time correlator Cf(t,t w ), which then gives the stationary cor- 
relator for t w — > 00. The detailed discussion will be presented 
elsewhere l49ll50ll . We start from the exact Eq. (O and use the 
projector a X y){a xy a xv y 1 {a xy as well as Eq. (O to get 



tw 



ds- 



xy 

) dC f (t,t w ) 



i m =0 

(30) 



<r(tw) = 7 / ds 
Jo 



{(7 X yG 0~ X y) 



(31) 



containing as numerator the familiar stationary shear stress, 
measured in 'flow curves' as function of shear rate IT371 l5a 
l57tl . A technical problem arises for hard spheres, where the in- 
stantaneous shear modulus (<r xy cr X y) diverges |3^l giving for- 
mally a = 0. The proper limit of increasin g st eepness in the 
repulsion has to be addressed in the future 14911 . In the spirit 

of the -model f58ll . we approximate the s-dependent nor- 
malized shear modulus by the transient correlator l34l[36ll . 



(&xye S(J xy) ^ g.f \ Geo 

\0 Xy O X y ] J 



(32) 



where we account for the different plateau heights of the re- 
spective normalized functions by setting Goo// ~ \', choos- 
ing a quadratic dependence $ 2 (s) would only change the re- 
sults imperceptibly. We will abbreviate a(t w — > 00) = a. 
The second factor on the right of Eq. ( f30b is the waiting 
time derivative, which we approximated in Eq. (l25b . We are 
hence able to show the two-time correlator for different wait- 
ing times for a glassy and a liquid state (Fig. [3J. The short 
time decay down onto the plateau / is independent of wait- 
ing time t w , whereas the long time decay becomes slightly 
faster with increasing waiting time. Overall the waiting-time 
dependence is small. In recent simulations of density fluc- 
tuations of soft spheres fiUl . the difference between the two 
correlators was found to be largest at intermediate times, and 
C{(t,t w ) < Cf '(f) was observed. Both properties are ful- 
filled by Eq. ( f30T >. Note that Eq. (f30b is exact in first order in 

tyj . 

Based on Fig. [3] and the knowledge about the transient cor- 
relators l36l l38ll . the short time decay of C{t,t w ) is inde- 
pendent of shearing for small shear rates. In glassy states at 
7^0 with yt — const., the transient correlator reaches 
a scaling function C[ • (7*) l36ll . and the two-time correlator 
from Eq. (130b reaches Cf(jt, 7^). 



VI. FINAL RESULTS FOR THE FDR 

Our final result for the susceptibility in terms of the waiting 
time derivative reads, 



Xf(f) 



(33) 



Eq. (1331 1 states the connection of two very different physical 
mechanisms: The violation of the equilibrium FDT and the 
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FIG. 3: C(t,t w ) from Eq. J30b for different waiting times as indi- 
cated by the different line styles. The model Eq. J28t is used 
to generate the transient correlators C^' (solid lines). We show a 
glassy state (e = 10 -3 ) at shear rates 7/r = 10 -8 (black) and 
j/T = 10~ 6 (red). Green and blue curves show the same respec- 
tive shear rates for a liquid state (e = — 1CP 3 ). The thin dashed line 
shows the equilibrium correlator C' e '(t) for the glassy state. For 
the liquid state, the equilibrium correlator coincides with the green 
curve. 



waiting time dependence of the two-time correlator at t w = 
0. The connection can be tested in simulations, where both 
quantities are accessible independently II 1 OL kill . The extra 
term in the FDT can indeed be connected to the time derivative 
of a correlation function reflecting its dissipative character, but 
no such simple relation occurs as in equilibrium. Using our 
approximation in Eq. ( fZ5b for the waiting time derivative, we 
can hence finally write our extended FDT 



Xl (t) « -C f (t) 




(') 



(*) 



Q (e) (t), 



(34) 



This equation connects the susceptibility to measurable quan- 
tities, at least in simulations, without adjustable parameter. 



A. FDR as function of time 



FIG. 4: C(t) from the i^^ -model with Eq. CP and x{t) via 
Eq. (34) for a glassy state (e = 1(T 3 ) and 7/r = 1CT 2 " with 
n — 1...4. Shown are integrated correlation, 1 — C(t) and re- 
sponse x (t) = J x(t')dt' ■ Inset shows additionally the transient 
correlator $ for comparison and the _X"( umv ) = I susceptibility for 
7/r = 10~ 8 . FromRef. O. 
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FIG. 5: Parametric plot of correlation C(t) versus response x'(t) = 
Jo x{t')dt' from Fig.|4](e = 10™ 3 ) together with constant non-trivial 
FDR (straight lines) at long times. The vertical solid line marks the 
plateau /. Inset shows the FDR X(t) as function of strain for the 
same susceptibilities. From Ref. (34ll . 



The relation between susceptibility and correlators requires 
correlators as input. We want to visualize the susceptibility 
using the schematic model d28l l. Eq. (f30b for the stationary 
correlator and Eq. d34b . Fig. [4] shows the resulting \ together 
with C for a glassy state at different shear rates. For short 
times, the equilibrium FDT is valid, while for long times the 
susceptibility is smaller than expected from the equilibrium 
FDT. This deviation is qualitatively similar for the different 
shear rates. For the smallest shear rate, we also plot x calcu- 
lated by Eq. (134-b with Cf replaced by Cf , see inset of Fig. |4] 
In this approximation, the FDR intriguingly takes the univer- 
sal value X( umv '(7t) = |, without any free parameters. The 
realistic susceptibility is achieved by including the difference 



between C^\t) and C(t). The parameter a is directly pro- 
portional to this difference. 

In the parametric plot (Fig. |5), the X = ^-approximation 
leads to two perfect lines with slopes — 1 and — | connected 
by a sharp kink at the non ergodicity parameter /. For the 
realistic curves, this kink is smoothed out, but the long time 
part is still well described by a straight line, i.e., the FDR is 
still almost constant during the final relaxation process. We 
predict a non-trivial time-independent FDR Xt ( jt) =const. if 

C f (and with Eq. (f30b also Cf) decays exponentially for long 
times, because A X f then decays exponentially with the same 
exponent. The slope of the long time line becomes smaller 
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FIG. 6: C{t) and x(t) via Eq. (SJ for a fluid state (e = -10~ 3 ) and 
■y/F — 10 _2n with n = 1...4. Shown are integrated correlation, 1 — 
C(i) and response x'(t) = J * x{t')dt'. Inset shows the parametric 
plot for the different shear rates. 

with increasing a (i.e., also with increasing value of Goo// 
in Eq. d32b). We find that the value of the long time FDR is 
always smaller than i in glassy states. 

The line cuts the FDT line below f for 7 — > 0. All these 
findings are in excellent agreement with the data in Ref. fioll . 
The FDR itself is of interest also, as function of time (inset 
of Fig. |5j. A rather sharp transition from 1 to | is observed 
when $ « C is approximated, which already takes place at 
jt m 10~ 3 , a time when the FDT violation is still invisible 
in Fig. [4] For the realistic curves, this transition happens two 
decades later. Strikingly, the huge difference is not apparent 
in the parametric plot, which we consider a serious drawback 
of this representation. 

Fig.|6]shows C{t) and \(t) f° r a fluid state. For large shear 
rates, these curves are similar to the glassy case, while for 
7^0, the equilibrium FDT holds for all times. In the para- 
metric plot (inset of Fig. |6]l one sees that the long time FDR 
is still approximately constant in time for the case 7 r _1 
(n = 3 in Fig- [6j, where shear relaxation and structural relax- 
ation compete, t is the a relaxation time of the un-sheared 
fluid. 

Summarizing, we find that the two separated relaxation 
steps lHH HH] (Fig. [3]) of the correlator in the limit of small 
shear rates for glassy states are connected to two different val- 
ues of the FDR. During the shear independent relaxation onto 
the plateau of height given by the non-ergodicity parameter /f , 

we have C f (t) (t) = C f (e) (t) + O^t) iU] in Eq. (O, and the 
equilibrium FDT holds. During the shear-induced final relax- 
ation from /f down to zero, i.e., for 7^0, and t — > 00 with 
t A f = const., the correlator without shear stays on the plateau 
and its derivative is negligible. A non-trivial FDR follows. In 
the glass holds 

r -ay) 7*<i, 

lim X ((t) = { -Cf(t) + it = 0(1), (35) 



FIG. 7: Long time FDR as function of shear rate for glassy (e = 
10"") and liquid (e = -10"™) states with n = 2, 3, 4. X(t -> 00) 
is determined from fits to the parametric plot as shown in Figs.f4]and 
[6] Inset shows the long time FDR as function of wavevector q for 
coherent (solid line) and incoherent (dashed line) density fluctuations 
at the critical density (e = 0). 

If one approximates stationary and transient correlator to be 
equal H, C f (t) (i) « C f (t), we find the interesting universal 
X = |-law for long times, 

Jim Yf(t c») = ~C f (t). (36) 
7^0 1 

The FDR, in this case, takes the universal value 
lim^o^fO -> 00) = X( univ )( 7 i) = |, independent 
of /. This is in good agreement with the findings in Ref. fioll . 
and corresponds to an effective temperature of T c ^/T = 2 for 
all observables. The initially additive correction in Eq. < fT3T > 
hence turns then into a multiplicative one, which does not 
depend on rescaled time during the complete final relaxation 
process. As summarized in Sec. U many spin models yield 
X = I at the critical temperature. The deviation from the 
value i of the long time FDR in our approach comes from 
the difference between stationary and transient correlators. 



B. FDR as Function of Shear Rate 

Fig. [7] shows the long time FDR as a function of shear 
rate for different densities above and below the glass transi- 
tion. The FDR was determined via fits to the parametric plot 
in the interval C(t) £ [0 : 0.1]. In the glass X(t — > 00) 
is nonanalytic while it goes to unity in the fluid as 7 — > 
(compare Fig. [6]). We verified that the FDT-violation starts 
quadratic in 7 in the fluid, as is to be expected due to sym- 
metries. X(t — > 00) is also nonanalytic as function of e and 
jumps to an finite value less than one. For all densities, the 
FDR decreases with shear rate. For constant shear rate, it de- 
creases with the density. This is also in agreement with the 
simulations fioll . 
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C. FDR as Function of Wavevector 

The realistic version of the extended FDT, taking into ac- 
count the difference of transient and stationary correlator, 
gives an observable dependent FDR in general. This can be 
quantified by using the exponential approximation for the long 
time transient correlator for glassy states (compare Eq. j29l ) 

C ( f t] (t -> oo) ps ft e -«fl-y|*. (37) 
The long time FDR then follows with Eqs. ([34]), ([30} and (f25b . 

Xf(t^co) = \ ~° f * . (38) 
1 — afcr 

The inset of Fig.[7]shows the long time FDR for coherent and 
incoherent density fluctuations at the critical density. We used 
the isotropic long time approximations (l29l and (1B241 > for af 
respectively and c a = 0.1 from Eq. ( [32b . The incoherent case 
was most extensively studied in Ref. lUOIl . The FDR in Fig. [7] 
is isotropic in the plane perpendicular to the shear direction 
but not independent of wave vector q, contradicting the idea 
of an effective temperature as proposed in Refs. [l^l and 
others. 

For q — > oo, h q j f q (corresponding to af) grows without 
bound and the FDR in Eq. ( [381 becomes negative eventually. 
For the parameters we used, the root is at q ~ 30. For larger 
values of ca (i.e., for larger values of Goo// in (1321)), the root 
is at smaller values of q. According to our considerations in 
the discussion section and the available simulation data, a neg- 
ative FDR is unphysical. For large values of q, the exponential 
approximation for the transient correlator or our approxima- 
tion (f3Qb for the two-time correlator might not be justified. 




10" 1 10° 10 1 10 2 10 3 
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FIG. 8: Comparison to simulation data for incoherent density fluctu- 
ations in the neutral direction (wave vector q = 7.47e z ) at tempera- 
ture T = 0.3 (T c = 0.435) and 7 = 10~ 3 . Circles and squares are 
the data (including units) from Fig. 1 1 in Ref. fich . lines are 1 — C q 
from Fig. 8 in Ref. fich . and the response x q W — Jo Xn(t')dt' cal- 
culated via Eq. 134t , The dashed line shows x q with approximation 

Cq'' ~ C q . Inset shows the different correlators, see main text. 

VII. HERMITIAN PART OF THE SMOLUCHOWSKI 
OPERATOR AND COMOVING FRAME 

In this section, we want to understand the violation of the 
FDT from a different point of view, i.e., from exact reformu- 
lations of the starting point (1131) . First, we split the SO into 
Hermitian and anti-Hermitian part to see how the Hermitian 
part is connected to the susceptibility in Eq. ( fT3b . We will 
then see that one can reformulate ( [T3l in terms of an advected 
derivative. 



D. Direct Comparison to Simulation Data 

Despite the dependence of the long time FDR on wavevec- 
tor, Eq. d34l ) is not in contradiction to the data in Ref. 1100, as 
can be seen by direct comparison to their Fig. 1 1 . For this, we 
need the quiescent as well as the transient correlator as input. 
Cq^ 1 has been measured in Ref. J59fl, suggesting that it can be 
approximated by a straight line beginning on the plateau f q of 
Cqit). In Fig. [8] we show the resulting susceptibilities. There 

is no adjustable parameter, when Cq « C q is taken. For the 

other curve, we calculated C q t} (t) by inversion of Eq. (130b . 
We used the dimensionless number a = 0.01 as fit param- 
eter which was chosen such that the resulting susceptibility 
fits best with the simulation data. The achieved agreement to 
X from the simulations is striking. In the inset we show the 
original C q from Ref. fioll together with our construction of 

C q e) and the calculated C q 4) . It which appears very reason- 
able compared with recent simulation data on Cf(t,t w ) J4lll 
and compared with Fig. [3] The value for f q used to construct 

Cq is also indicated in the main figure. 



A. Hermitian part 

Investigating the stationary correlator in Eq. (0, one finds 
that the operator Sv is not Hermitian in the average with $> s 
i60ll . This is why one cannot show that Cf (t) is of positive 
type l6lll via, e.g. an expansion in eigenfunctions since only 
an expansion into a biorthogonal set is possible JHH^]. 

Subsequent to realizing this, we want to split the SO into 
its Hermitian and its anti-Hermitian part with respect to the 
average with $> s . Recall that is the adjoint of S¥ in the 
unweighted scalar product U, H3, 0]. The adjoint of fit in 
the stationary average is defined by 

(sfi f f} w = J dr*.(r) ff (r)fit/*(r) 
= J aT* s (r)/*(r)f>t.9(r) = (/*fit 5 ) 7 . (39) 

We already stressed that fit is neither identical to Q, nor 
ft = J2 9 i+ 2 ^ ln *•) • °i - F < ' d i -di K Ti. (40) 

i 
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The difference between non-equilibrium forces di In '5 s and 
the potential forces Fj = diln^ e appears |{48tl . Now the 
Hermitian and the anti-Hermitian parts of fit with respect to 
stationary averaging are given by 

fit, =i(fit + fit) = df + (di In*.) • ^ , (41a) 

i 

fi^=i(fit-fit) 

= ^ -{di In * s ) • 9, + Fi ■ d t +8i-K- Vi . (41b) 

i 

We obviously have fit = £V H + SV A . SV H is similar to the equi- 
librium SO fit with forces F j replaced by the non-equilibrium 
forces di In * s . As expected, the anti-Hermitian part contains 
the shear part <5fit. It also contains the difference between 
equilibrium and non-equilibrium forces. The eigenvalues of 
Q A are imaginary and the eigenvalues of fij, are real JH]. In 
the given case, fij, can furthermore be shown to have nega- 
tive semi-definite spectrum as does the equilibrium operator, 
because we have 



We hence finally have 



(7) 



dr l dr l / 



(42) 



If the correlation function is real for all times, Cf(t) = C f * (£), 
as can be shown e.g. for density fluctuations |38|. the initial 
decay rate is negative since il A does not contribute, 



K 



(7) 



(7) 



(43) 



Thus a real correlator initially always decays, i.e., the external 
shear cannot enhance the fluctuations. Higher order terms in 
t contain contributions of fi^ and such an argument is not 
possible. 



B. Susceptibility and comoving frame 

We now come to the connection of fi^ to the susceptibility. 
Eq. dT3b can be written 



Xfg 



(7) 



(44) 



The response of the system is not given by the time derivative 
with respect to the full dynamics but by the time derivative 
with respect to the Hermitian, i.e., the "well behaved" dynam- 
ics. It follows that we can write the susceptibility, 



X is = -(r(rt-n A )e nU g) 



(7) 



(45) 



We note that this equation is very similar to Eq. (13) in 
Ref. H62J1 since h(L — L*) is the anti-Hermitian part of L in 
Ref. 116211 . Eq. (l44l can be made more illustrative by realizing 
that fi^ can be expressed by the probability current jf, 



fi T 



(46) 



Xfg 



(7) 



(47) 



The derivative in the brackets can be identified as the con- 
vective or comoving derivative which is often used in fluid 
dynamics J63ll . It measures the change of the function in the 
frame comoving with the probability current. If one could 
measure the fluctuations in this comoving frame, these would 
be connected to the corresponding susceptibility by the equi- 
librium FDT This was also found for the velocity fluctuations 
of a single driven particle in Ref. 0], The difference in our 
system is that the probability current, i.e., the local mean ve- 
locity, speaking with the authors of Ref. |7|], does not depend 
on spatial position x, but on the relative position of all the 
particles because it originates from particle interactions. 

Let us finish with interpreting the comoving frame. Xf de- 
scribes the tendency of particles to move with the stationary 
current. If the stationary current vanishes, we have Xf = 1. 
If the particle trajectories are completely constraint to follow 
the current, we have Xf = 0, because a small external force 
cannot change these trajectories and \f = 0. As examples 
for the latter case, let us speculate about the experiments in 
Refs. ||64L l65tl . They consider a rather dilute suspension of 
colloids in a highly viscous solvent. The bare diffusion coef- 
ficient is approximately zero (so called non-Brownian parti- 
cles), i.e., on the experimental timescale the particles do not 
move at all without shear. Under shear, the particles move 
with the flow and one observes diffusion in the directions per- 
pendicular to the shearing due to interactions. A very small 
external force does not change the trajectories of the parti- 
cles (on the timescale of the experiment) due to the high vis- 
cosity. We expect Xf = in this case because the parti- 
cles completely follow the probability current. The studies 
in Ref. 11641 16511 do not consider the susceptibility, the focus 
is put on the question whether the system is chaotic or not. 
The finding that the dynamics is irreversible under some con- 
ditions makes it even harder to predict the FDR, which would 
be of great interest. 



C. FDT for eigenfunctions 

From Eq. (O and fi^ = |(fit + fit), we find for arbitrary 

/ = f({xi,Vi,Zi}), 

»(*)=-|(m t e nt «/) W) -|((n t ne nt ^ W) 



(48) 



(fiWV 



This form is especially illustrative since it explicitely shows 
that the FDT violation occurs because fit is not Hermitian in 
the stationary average. If it was, the two terms above would be 
equal and the equilibrium FDT would hold. We note that this 
form is equivalent to Eq. (1 1) in Ref. rf62tl . As pointed out by 
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Baiesi, Maes and Wynants, Eq. (f48b in the case of simulations 
has the advantage that correlation functions of well defined 
quantities (/ and fit /) can be evaluated. This indicates the 
usefulness of Eq. d48l l relative to Eqs. ( TTBT l and d47b . 

If we consider the case that / = <j> n with tp n eigenfunction 
of fit, fit0 n = X n (f> n , we find 



The equilibrium FDT thus holds for f = q 



(49) 




FIG. 9: The shear step model: The particle is trapped in wells of 
width a with infinitely high potential barriers. The center-to-center 
distance of the wells is denoted b. 



VIII. DISCUSSION 



A. Deterministic versus Stochastic Motion 



We saw in Sec. IVIII that the susceptibility measures the 
fluctuations of the particles in the frame comoving with the 
probability current jf . We conclude that we can split the dis- 
placements of the particles into two meaningful parts. First, 
the stochastic motion in the frame comoving with the average 
probability current. Second, the motion following the average 
probability current, which is deterministic and comes from the 
particle interactions. The deterministic part is not measured 
by the susceptibility, \ is thus smaller than expected from the 
equilibrium FDT. It measures only parts of the dynamics. We 
have X{ < 1, Let us quantify the above discussion as far as 
possible. In Eq. d44b . we see that we can formally split the 
time derivative of the stationary correlation function into two 
pieces, the stochastic one, measured by the susceptibility, and 
the deterministic one following the probability current, 







(7) 



me 



(7) 



The inequality Xf < | in glasses, see Fig.|7]translates into an 
inequality for the two derivatives above, 



(-) 



> 



f fike n \ 



(7) 



(50) 



In completely shear governed decay of glassy states, the deter- 
ministic displacements of the particles due to the probability 
current are larger than the stochastic fluctuations around this 
average current. In other words, if the stochastic motion was 
faster than the deterministic one, the decay would not be com- 
pletely shear governed. 

It is likely that with increasing density or lowering temper- 
ature, the particles are more and more confined to follow the 
probability current and the FDR gets smaller and smaller and 
might eventually reach zero. 



B. Shear Step Model 



simplify the dynamics of super-cooled liquids and glasses, be- 
cause the particles themselves form the traps for each other 
and it is therefore not easily possible to map the problem onto 
a single particle problem. Nevertheless, we want to introduce 
a simple toy model which will provide more insight into the 
FDT violation of sheared colloidal glasses. The model is de- 
picted in Fig. [9] The particle, surrounded by solvent (diffusiv- 
ity Do = ksTfiQ, we restore physical units), is trapped in an 
infinite potential well of width a. This potential well is one 
in a row of infinitely many with center-to-center distances b. 
The particle diffuses in the well and a very long time after we 
measured its position, its probability distribution P(z) is con- 
stant within the well. The shear, i.e., a mechanism which lifts 
the particle over the potential barriers is introduced as follows; 
At time t = 0, 20, 30, . . . , the 'shear steps' lift the particle 
into a neighboring well according to its current position: If it 
is on the left side of the well, it gets to the well on the left 
hand side, if it is on the right hand side, it gets into the well 
on the right hand side. The model hence describes a direction 
perpendicular to the shear flow, e.g. the z-direction, where 
the influence of shear is symmetric. The initial position of the 
particle in the new well shall be distributed randomly. The 
resulting dynamics has some similarities to a Cauchy process 
lf72tl and shares many properties of colloidal suspensions as 
will be shown below. 

Detailed balance is broken because the reverse step, that a 
particle is taken out of the right side of a trap and put back 
into its left neighbor, is missing. 

The shear step model is much simpler than other trap mod- 
els considered in the literature. The Bouchaud trap model Irolll 
contains a distribution of traps of different depth, allowing to 
study different situations such as aging. The simplicity of our 
model makes it easier to be analyzed and the result for the 
FDR contains only the parameters a and b, whose values are 
of comparable size. 

We regard the limit of small shear rates, it corresponds to 
^> a 2 /Do, i.e., the time between two shear steps is much 
longer than it takes the particle to relax in the well. We will 
first present the mean squared displacement of the particle and 
then its mobility under a small test force. 

At short times, the mean squared displacement (MSD) of 
the particle is the one of a free particle, 



Trap models have repeatedly been used to study the slow 
dynamics of glassy systems and to investigate the violation of 
the equilibrium FDT H M, HI M, EZS Hfl] . They boldly 



\im((z - z(O)Y) =2D t. 
t— >o 



(51) 



For times > t ^> a 2 /Do, the dynamics of the particle is 
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FIG. 10: Mean squared displacement without (red squares) and mean 
traveled distance with external force F (blue circles) in the shear step 
model with b = a and O = 2a 2 / (2Do). Solid line shows the short 
time asymptote of Eqs. i5U and i54\ . Dashed line shows the plateau 
value of ^p, Eqs. l !52t and J57b , The dotted line shows the long 
time asymptote from Eq. d53l >. dash-dotted line the asymptote from 
Eq. (39). 



glassy, i.e., the MSD is constant on the plateau. This plateau 
value can be derived from the constant probability distribution 
of the particle in the well, P(z) = 1/a, 



((z-z(0)f 



1 



a/2 



a/2 



dz(0) 



a/2 



a/2 



dz(z-z(0)) 2 = 



6 ' 
(52) 

We notice that in accordance with the glassy dynamics of col- 
loidal suspensions, the plateau value is independent of Dq and 
temperature. Note that the initial positions are distributed with 
P(z) as well. At long times, the particle performs a random 
walk with step-length b and number of steps t/0 lf73ll and the 
MSD approaches 



((z-z(0)Y 



,t_ 

0' 



t»e. 



(53) 



The long time dynamics is independent of temperature and 
Dq, as is the long time decay of the density correlator for 
sheared colloidal glasses l36l 13^1 . The timescale is set by 
0, corresponding to the a 'shear rate' of 7 = _1 . 

To calculate the fluctuation dissipation ratio X, we have 
to find the mobility fi(t) in response to a small test force F 
acting on the particle, say to the right. This test force shall 
not influence the jump rules defined above. At short times the 
particle obeys the (integrated) equilibrium FDT (the Einstein 
relation), 



lim — ((z 

t^Q F u 



'(0))> 



1 



lim ■ 
t^o 2k B T 



((z - z(0)) 2 ) = fi Q t. 

(54) 



(. . .y' denotes an average under the influence of the exter- 
nal force. The external force changes the probability distri- 
bution of the particle in the well. It is more likely to find the 



particle on the right hand side of the well than on the left hand 
side. The distribution P F (z) for 9 > t 3> a 2 / D follows 
the Boltzmann distribution i43ll . which, in linear order of the 
external force, reads, 



P F (z) 



[ h(i + T5t) -a/2<z<a/2, 
\ else. 



(55) 



The mean traveled distance is easily derived, 



i{(z-z(0))) (F) 



0/2 



a/2 



dz(0) 



a/2 



a/2 



dzP F {z)(z - z(0)) 



Ylk B T' 



(56) 
(57) 



Comparing with Eq. ((52), the equilibrium FDT holds for all 
times t < as expected. 

Due to the distorted probability distribution Pf{z) in 
Eq. dB3T l, the shear step at t = Q will take the particle more 
likely to the right. The rate R for jumps to the right minus the 
rate for jumps to the left follows, 



R = 



J a/2 P F {z)dz - J°_ a/2 P F {z)dz 

rJ?Mz)dz 



4 knT 



0(F 2 ). 
(58) 

At every shear step, the particle travels on average the distance 
bR. For t > 6, the initially traveled distance in the well is 
negligible and we have 



i{(z-z(0))) (F) 



a b t 
4k B TO' 



t»e. 



(59) 



The mobility of the particle is finite at long times and inde- 
pendent of the diffusivity Do. The FDR X is in this situation 
defined by 



((z-z(0)) 2 ). (60) 



We illustrate the results above with simulations of the de- 
scribed dynamics, see Appendix |Cl for details. Fig. [TOl shows 
the simulation results together with the derived asymptotic 
formulae. We see that the FDT holds for times t < <3 and 
is violated for t > 8. 



1. Long Time FDR 

We see in Fig.[l0]that the Einstein relation does not hold for 
long times. The long time FDR X(t 0) is different from 
unity and follows with Eqs. (|53l and 1 



X(t) 



1 t<e, 
-- t ^> 

2 b L ^ ^1 



(61) 



This is illustrated in Fig.fTTl 
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FIG. 1 1 : Parametric plot of MSD versus mean traveled distance for 
b = a and O = 2a 2 /(2Do). Solid line shows the short time relation 



X — 1. Dashed line shows the long time FDR X = 
the same graph zoomed into the short time part. 



|, Inset shows 



2. Discussion 

Let us emphasize the similarities between the shear step 
model and the sheared colloidal glass. For both, X(t 6) 
is independent of temperature because the potentials are in- 
finitely high. In soft sphere glasses, this is not true because 
the temperature also governs the glass transition. 

In the colloidal glass, the particle is trapped in the cage. The 
short time motion, the rattling in the cage, is FDT-like, as is 
the motion of the particle in the potential well. For long times, 
the shear drives the particle into a neighboring cage. For the 
directions perpendicular to the shear, this motion is symmetric 
(without test force) as in the shear step model. In the shear 
step model, the small external force does not change the shear 
step mechanism, it only changes the probability distribution 
of the particle in the well and thereby the motion becomes 
asymmetric. It is appealing to imagine something similar to 
happen in the colloidal system: The external force influences 
the distribution in the cage and makes the probability for the 
particle to be driven to the neighboring cages asymmetric. In 
both cases, this mechanism leads to a finite mobility which 
does not obey the equilibrium FDT 

In the shear step model, X(t > 9) reaches | if b goes to 
a. b < a (leading to 1 > X > |) is physically not reasonable 
and we get the constraint of X < | which coincides with the 
one for the real system, see Fig. [7] 

The model allows to discuss two more effects: In the col- 
loidal system at increasing density, the cages become smaller 
and smaller. This corresponds to a decreasing value of and 
X(t 3> Q) decreases. This is in accordance with Fig.|7]and 
Ref. flQl . In the shear step model, it eventually reaches zero 
being always positive. 

For increasing shear rates, i.e., decreasing 0, the particle 
has less time to relax in the well and the effect of the external 
force decreases. Fig. [T2l shows that decreasing 6 lowers the 



FIG. 12: Long time FDR in the shear step model as function of 'shear 
rate' 1/6. 



value of X. It is also in agreement with Fig.|7]and the simula- 
tions in Ref. fioll . We believe that the decrease in the colloidal 
system has the same origin, i.e., the particle has less time to 
adjust its distribution in the cage in response to the force. 

The FDT has been studied in Bouchaud's trap model pre- 
viously ||7b"ll . where in contrast to our model, the parametric 
plot has continously varying slope. The model in Ref. 17011 is 
more realistic than ours and different observables can be stud- 
ied. The advantage of the shear step model is its simplicity 
and its time independent long time FDR which has only the 
physically illustrative free parameters a and b. 



IX. SUMMARY 

We investigated the relation between susceptibility and cor- 
relation functions for colloidal suspensions at the glass transi- 
tion. While the equilibrium FDT holds at short times, a time- 
independent positive FDR smaller than unity is obtained at 
long times during the shear driven decay. We find that the 
long time FDR is nearly isotropic in the plane perpendicular 
to the shear flow and takes the universal value X{-yt) = i 
in glasses at small shear rates in the simplest approximation. 
This agrees with the interpretation of an effective temperature. 
Nevertheless, corrections arise from the difference of the sta- 
tionary to the transient correlator and depend on the consid- 
ered observable. They alter Xf to values Xf < \ in the glass. 
Our findings are in good agreement with the simulations in 
Ref. Hi. 

While we used as central approximation the novel relation 
for the waiting time derivative, Eq. d25b . a more standard MCT 
and projection operator approach leads to qualitatively equiva- 
lent results; see AppendixlAl Considering the crudeness of the 
MCT decoupling and of our approximation (l25l l. the quanti- 
tative differences between the two approaches appear reason- 
able. From both approaches we can conclude that there is a 
nontrivial FDR Xf (7*) = Xf for small shear rates during the 
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final shear driven decay. 

X[ depends on shear rate non-analytically for all shear- 
melted glassy states. At the glass transition, the 7^0 value 
jumps discontinuously from its nontrivial value X{ (7 — > 0) < 
1 to the equilibrium value Xf = 1. For finite drive, Xf 
decreases below unity for all states. The discontinuous be- 
havior of Xf (7 — > 0) results from the shear driven decay on 
timescale 7 -1 . Within MCT-ITT the shear governed final de- 
cay is also the origin of a finite dynamic yield stress which 
also jumps discontinuously to its equilibrium value at the glass 
transition. These predictions differ from the mean-field spin- 
glass results. Ref. ifioll finds power-law-fiuid behavior but no 
dynamic yield stress. Moreover, the FDR at vanishing shear 
rate moves continously to the equilibrium value unity at the 
glass transition. The MCT-ITT scenario for yielding and fluc- 
tuation dissipation relations is thus unique compared to other 
approaches to shear-melted glasses. Investigation of shear 
driven glassy dispersions thus provides a unique possibility to 
discriminate between different theories of the glass transition. 

The incoherent motion and the relation between diffusivity 
and mobility can also be studied within our approach, and is 
topic of a companion paper lp74tl - 

Our finding of values close to the universal X(jt) = ^ 
points to intriguing connections to critical spin models. Open 
questions concern establishing such a connection and to ad- 
dress the concept of an effective temperature, which was de- 
veloped for ageing and driven mean field models. It might 
also be interesting in the future to study the response of the 
system to small perturbations in the shear rate 7 17511 . 

We derived a relation between the dominant part of the vio- 
lating term and the waiting time derivative of Cf (t, t w ), viz a 
relation of two completely different physical quantities. This 
connection can be tested in simulations. Because we identi- 
fied all but one contributions of the violating term with inde- 
pendently measurable quantities, all our approximations are 
testable independently in simulations. The difference of the 
measured A\i in FDT simulations to the sum of the measured 
terms of A\i in waiting time simulations yields the contribu- 
tion of the last term in A%f. 

We presented a new exact formulation of the susceptibility 
which involves the Hermitian part of the SO. This part is inter- 
preted to represent the dynamics in the frame comoving with 
the probability current. 

We introduced the shear step model to illustrate the FDT 
violation in sheared colloidal suspension. In the shear step 
model, the FDR takes values X < 5. 
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APPENDIX A: MODE COUPLING APPROACH 

In this Appendix, we start over and present the analysis 
of the susceptibility in Eq. ( fT6l l for density fluctuations using 
Zwanzig Mori projections, f = Qq = J2i e* q r< for coherent 
and / = £>q = e iqrs for incoherent density fluctuations. We 
will treat the two cases at once with g denoting either £ q or 
£q. / being x-independent translates into q x = 0. Normalized 
equilibrium, transient and stationary correlators, as defined in 
SecEAl are denoted $^ e) (i) = (g* e n - 1 g) / (g* g) , $ q (f) = 
(g*e nH g)/(g*g) and G q (t) = (g* e nU g>« / (g* q)W . Sta- 
tionary averages are normalized with (g*g)^\ the initial 
value of the stationary correlator. In the coherent case, this 
is the distorted static structure factor Sq = (f? q £q)^ /N, in 
the incoherent case, it equals unity, g q *£> q = 1. The transient 
coherent correlator is normalized by the equilibrium structure 
factor S q = (g*Q g )/N. The normalized violating term is then 
defined by (compare Eq. $13[ ) 

a ^ (^H (^(?£■ 8 ' en, ' !, ) tt, + <; '' (^, ■ 

(Al) 

One gets the normalized analog to Eq. ( IT61 ). 

i ' 

Since we are left with an equilibrium average, it is useful to 
express Axq{t) in terms of the transient correlator $ q (i) as is 
done in the following subsection. 

1. Zwanzig-Mori Formalism - FDT Holds at t = 

We use an identity obtained in the Zwanzig-Mori projection 
operator formalism ll76tl (see Eq. (1 1) in Ref. ll77lo to find the 
exact relation 

A X q(t) = / dt'Nq(t-t')$q(t'), (A3) 

Jo 

Q e Q^QtQ[H e \ (A4) 

with P = g)(g* g)^ 1 (g* projecting on a subspace of density 
fluctuations and Q = 1 — P. Nq(t) can also be split into the 
three contributions according to Eq. < fT~8T >, which will be done 
at the very end only. Eq. dA3b could be expected to contain a 
second term, A\q(t = 0)<I> q (t), the static coupling at t = 0. 
It vanishes in (TAli i.e., Ax q (0) = 0; Ax q is real JH], and 
second and third term in Eq. < fT~8T > cancel at t — 0. The waiting 
time derivative vanishes at t = due to symmetry for arbitrary 
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(<r xy f*f) = 0. 



(A5) 



It has thus been shown that the equilibrium FDT is exactly 
valid at t = 0. 



2. Second Projection Step 

MCT approximations for the function AT q (i) directly are 
not useful because they cannot account for the fact that Nq(i) 
is a fast function. This is achieved with a second projection 
step following Cichocki and Hess lf78"ll . see also Refs. 0411361 
HH. The adjoint of the Smoluchowski operator is formally 
decomposed as 



(A6) 



The function AT q (t) is then connected to a new function n q (t) 
governed by fi*. The following identity can be proven by dif- 
ferentiation, 

JV q (t)/r, = «,(*)- / dt'Nq{t')mq(t-t'), (A7) 
Jo 



9r, 



m q (i) 



(A8) 



(A9) 



T g = — (q*SVq)/(q*q) is the initial decay rate. It equals 
q 2 1 ' S q for the coherent and q 2 for the incoherent case. m q is 
identified as the memory function, which appears in the equa- 
tion of motion of the transient correlator. It is related to the 
transient correlator for q x = exactly by l37ll 

$ q (t)+r„ j* q (t) + J dt'm^it - 0*q(«')} = °- (A10) 

We will later be interested in the Laplace transform, m q (z) = 
/ °° dte~ zl ! m q (i), which at z = is in the limit of slow dy- 
namics given with Eq. ( lAlOt by lf79ll . 



(All) 



m q (z = 0) = $ q (z = 0). 



The benefit of the second projection step can now be illumi- 
nated by regarding Nq(z), which is given via Eq. ( |A7t by 



Nq(z) = 



n q (z) 



l/T q + rriq(z)' 



(A12) 



To discuss this, we turn to scalin g co nsiderations. m q (i) de- 
cays on timescale (t _1 + 7|) _1 113611 . where r is the a relax- 
ation time of the un-sheared system. In the glass, it is formally 
infinite and we have 



lim m„(z = 0) = 

7— fO 



0(t) liquid, 
O(^) glass. 



(A13) 



As will be shown below (Eq. (IB 171 )) by MCT approximations, 
Tiq(z) has the following properties, 



lim n q (z = 0) 
7^0 



0(7 2 r 2 ) liquid, 
O(l) glass. 



(A14) 



The scaling ?n q (z = 0) oc I7I 1 in the glass provides Nq(z = 
0) via Eq. ( IA12I ) with an additional power of 7 compared to 

n q (z = 0), 



lim NJz = 0) 

7^0 



C(7 2 r) liquid, 
C(|Tl) glass. 



(A15) 



This is the benefit of the second projection. Note that N^{z = 
0) corresponds to a rate, which is always small but changes at 
the glass transition. Because of its smallness, Nq(z = 0) is 
difficult to approximate quantitatively. 

The time integrated violating term Ax q finally follows at 
small shear rates with Eq. 



liquid, 
glass. 



(A16) 



Eq. (IAI6I 1 is in accordance with simulations and with the 
physically expected property of x q to be always finite at 
z = 0. The response should not diverge. 

In both the liquid and the glass, the violating term Ax q is 
symmetric in 7, reflecting the fact that fluctuations in z- and in- 
direction are independent of the direction of shearing. While 
Ax q is analytic in 7 in the fluid, it is nonanalytic in the glass. 

We can now compare the derived property in Eq. ( IA16b to 
the approximation in Eq. d34l i. from which we find for glassy 
states and 7^0, 



A Xf (z = 0) 




Cf(t, t w ) 



t w =o 



(* = 0) « -\f u 



(A17) 

in accordance with Eq. (IA16b - fi is the non-ergodicity param- 
eter. In the fluid, Eqs. (IA16b and d34b also agree which can 
only be shown numerically, see Fig. [7] 



3. Markov Approximation - Long Time FDR 

Using a special projection step, we have shown in the pre- 
vious subsection that the function Nq(z = 0) is of order I7I in 
the glass, i.e., we have reason to assume that A r q (<) decays fast 
in time compared to < E > q (t) which diverges like | -^y- 1 1 at z = 0. 
With this assumption, Eq. dA3b can be written in Markov ap- 
proximation using the (5-function, iV q (t) w iV q (z = Q)S(t). 
For the susceptibility follows 



Xq(*) 



d_ 

dt 



Cq(t) + Nq( Z = 0)$ q (t). 



(A18) 



We will see in Subsec. lA4 al that Eq. ( IA1 8b gives very similar 
results to Eq. (f34-b . According to Eq. (IA18b . the equilibrium 
FDT is violated if Nq(z = 0) is nonzero. For short times, 
jt < 1, we have \Nq(z = 0)[ < \Cq(t)\ and the equilib- 
rium FDT holds. Nq(z = 0) is of order |-y [, see (IAT51) . Cq(t) 
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is of order T q = 0(1). For long times, jt w 1, C q (t) is 
also of order the two terms are comparable in size and the 
equilibrium FDT is violated. The long time FDR is addition- 
ally independent of time, if C q (i) and $ q (i) decay exponen- 
tially for long times with the same timescale. Approximating 
the two correlators to be equal and relaxing exponentially in 
glassy states, 



lim C q (f) « $ q (t) M f q e 



-Oq|7|i 



we find in the glass, 



Xq(*) = 



(A19) 



(A20) 



This equation is in qualitative agreement with Eq. ( l35l l. The 
long time FDR 



XJt 



N^{z = 0) 



(A21) 



is time independent and also independent of shear rate for 
7^0. It is hence non-analytic as pointed out before. As 
we will see, N q (z = 0) is negative in MCT approximations 
and the FDR is smaller than unity in agreement with Eq. ( [34-b . 



4. FDT Violation Quantitative - Connection of the two 
approaches 

We will now perform MCT approximations for the memory 
function n q (i). It contains two evolution operators, one for 
the correlation time t and one for the transient time s, which 
entered through the ITT approach. We rewrite Eq. (IA81 > via 
the identity (JT7J which leads to three contributions for n q (£) 



(2) 



Mi 



d S (cr xy e n ' s ^g*U(t)^g), 



7 



(a xy g*U(t)^g) 



J ds(a xy e n1s 



wffl(t)= — ; / ds(<J T ,,e 



n 1 & 



U{t) =Qe Q(VQt Q. 



(A22a) 



(A22b) 



(A22c) 



The ,s-integration in n q ^ could be done directly as in Eq. (TT~9b . 

Also ra q 2) (i)and n q 3) (i) can now be identified with derivatives 
with respect to s. It is important to note that without identify- 
ing these s-derivatives, the correct 7-dependence of A% q (i) 
would not be achieved. 



The detailed MCT approximations for the terms above are 
shown and evaluated in Appendix [B] In the following sub- 
section, we compare the Appendix approach to the main-text 
approach. 



a. Connection of the two approaches 

Let us compare the three terms of A\i in Eq. ( TT~8T > to 
Eq. dA22b . From the analysis in this Appendix, we have the 
exact relation for f = g, which is more handy in Laplace 
space, 



n q (z) 



1/T q + TO q (z) 



*q(*)- 



(A23) 



The waiting time derivative in Eq. ( fT8l ) is for density fluctua- 
tions exactly given by 



*™=o 



l/Fg + TO q (z) 

(A24) 



And the sum of the other two terms is exactly given by 

'£,(*)- *,(*)] + A Xq 3 >(t)}(z) 



4 



1/T q + m q (z) 



* q (z). (A25) 



With the Markov approximation in Eq. ( |A18t , these are sim- 
plified to (note that m q (z = 0) oc l^l" 1 > l/r g ) 



1 _d 

2 dt 



A 1) (z = 0) 
m q (z = 0) 



*q(*)> 



(A26a) 



{^23} (*) « " g)(z = 0) / + W !l (z = 0) ^(0. 



m q (z = 0) 



(A26b) 



with n q ^(z = 0)/m q (z = 0) oc |-y| as 7 — > 0, see Eq. ( IA14I ). 
From Eq. $25[ , which we consider very accurate, the wait- 
ing time derivative in glassy states is for long times equal to 
$ q (£). We conclude that Eqs. d25l l and ( |A26at are in qualita- 
tive agreement, if the transient correlator decays exponentially 
for long times with timescale oc |7| _1 . This holds well Oal 
and also has been used in the main text, see Eq. (|29| |. 

In Fig. [2] we show the quantitative comparison of the func- 
tions at z = 0, see Appendix [B] for details on the MCT ap- 
proximations for the terms ( lA22b and their numerical evalu- 
ation. We see, that the MCT-estimates compare quite well to 
the prediction of Eq. ( |25T l. Fig. [T"3l shows the long time FDR 
as function of q as calculated by the MCT approximations for 
?i q and Eq. ( IA2U . We show only the contribution of the first 



term, i.e., n q (z = 0) 



(z = 0). Also, in Eq. dA2fl . 



the difference between stationary and transient correlators is 
neglected since it follows with Eq. ( |A19b . According to our 
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z-direction inc. 
z-direction con. 
y-direction inc. 
y-direction con. 

15 

q 



30 



At the left hand side appears the expression 

y(l) _ (^jygqQgkgq-k) 

qk = NSg 

S' S' 
« k x {k v -q y )-^S k + k x k y ^-S q . k . (B4) 

On the right hand side, we have the standard vertex 13 



Vg = 



(gkgq-kQ^gq) 



NS 



k-qdk 



q ■ ((k - q)nck~ q + kriCfc). 



(B5) 

In the derivation of Eqs. dB4b and dB5l ). the convolution ap- 
proximation for the static three point correlation function was 
used JH, 



FIG. 13: Long time FDR as function of wavevector q for density 
fluctuations in z- and y-direction in the limit of small shear rates. 
The packing fraction is the critical one (e = 0). 



analysis in the main text, Eq. d25l l. these two simplifications 
would yield X = | for all q. We see that the FDR evalu- 
ated from dA21b depends rather strongly on wavevector, but 
given the complexity of the involved functions, the result is 
still satisfying. 



APPENDIX B: MCT-APPROXIMATIONS FOR n q (i) 

Here we show the detailed approximations for the formally 
exact expressions for n q in (IA22b by projection onto densi- 
ties. This physical approximation amounts to assuming that 
these are the only slow variables, sufficient to describe the re- 
laxation of the local structure in the glassy regime. 



(gqgk-qgk) ~ NS q S k -q S k 



(B6) 



(2) (3) 

The treatment of n q ' and rt q ' is more involved. We first sep- 
arate the time evolution operator in t with pair projectors to 
get, 



7 



2NS£ ) T q Jo 



J ds(a xy t 



Ql^QP 2 U{t)P 2 {t)Q9) eq ), (B7) 



< 3) (*) 



2NS$ ) T q Jo 



ds(a xy e n ^ 



(rtg*)QP 2 U(t)P 2 (t)Qrt eci 



(B8) 



After doing this, we are on the left hand side of the projectors 
left with the two respective expressions, 



(B9) 



1. Coherent Case 

For coherent density fluctuations, the time dependent pair 
density projector is given by J38ll 



m = E 



k>p 



The memory function n q has only one time evolution opera- 
tor and can hence be approximated via 'standard' routes with 
the projector P 2 {t) IBal . 



gk(t)g P (t))(gk(t)gp(t) 
N 2 S) 5 _(t)S p ( t ) 



(Bl) 



nW(t) 



7 



2NS^T q 



(B2) 

For the appearing time dependent four point correlation func- 
tion, the factorization approximation is used 



(gk'^gk'C-tJ-q^Wgkgk-q) « 
N 2 S k{ _ t )S k (-t)- q <S>k(-t)(t)<f>k(-t)-q(t)5 k ,V 



(B3) 



and 



(T 2; j / e 0ts (ri t g q ) Q £»q-kgk 



(BIO) 



Writing Q as 1 — P, we realize that the term containing P is 
identical for both terms (they are real), 



<7 X yC 



(^g q ) gq 



1 



NS q xeq 



(gqgq-kgk), (Bll) 



with opposite sign. These terms cancel each other. We are left 
with the two expressions, 



and 



<Tx V e Q s g q ^ t g q -kgk 



CT:zye nts (jVp q )£> q _kgk 



(B12) 



(B13) 



There is in principal more than one option to treat these terms, 
but we will argue that only one option is applicable. The stan- 
dard way, i.e., the usage of P 2 right and left of the time evo- 
lution operator is not preferable since it would not preserve 
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the derivative with respect to s. As already mentioned, this 
derivative is necessary for the correct 7-dependence. That is 
why we chose to use the triple densities projector P3, 



Ps(t) 



E 

k>p>n 



0k(t) Qp(t) Qa{t) ) (^(t) B P (t) 8n(t) 



(B14) 



Eq. (IB 1 21 > is written as 



Pq-kPk 



0- :E yf3(-s)e 0ts e*ilVq-k£'k 

(B15) 

We have to demand that the wave-vectors in the triple pro- 
jector take the values of the wave-vectors on the right hand 
side. Due to this constraint, the summation in Eq. ( IB 14b con- 
tains only one term and no counting factor appears. The left 
hand side is the Vertex in Eq. (IB4I ) for s-dependent 
wave-vectors (the projector Q does not make any difference in 
( lB4b ). The appearing six point s-dependent correlation func- 
tion is approximated as, 



(&J0k(_ a )0k(-«)-q e 



n 1 1 



£>k£>k-q) 



N 3 S q S k{ . 



-qS k (- 



*q(«) 



$k(-a)(»)$k(-»)-q(f 



This approximation rests on the observation that the operator 
fit acts as an s derivative. $j c (-s)( s ) an d ^k(-s)-q( s ) depend 
on s via the decay of the correlator and via the s-dependent 
wave-vectors. Since fit in ( IB 16b represents the s-derivative 
with respect to correlator dynamics, we have to subtract the 
change in s due to the change of the wave-vectors. The term 
in Eq. (IB 1 31 > is treated analogously with the triple density 
projector. Here, the approximation for the appearing six point 
correlation function is more straight forward, since $ q (s) has 
no wavevector advection, 
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(fitpq^kgk-q) 



N 3 S a S, 



q £>k(—s)—q £>k( — s) 
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^$q(s)J $k(- S )0)$k(- 8 )-q(s)- 



Collecting the terms, we finally find the following expres- 



sions, where n^* 1 and n^' denote the functions without the 
terms in Eq. (IB 1 lb . 
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(B16a) 

$k(-.) (»)*k(-«)-q(*)) $k(t)*k-q(*). 

(B16b) 



(B16c) 



n—k 



rt "X rL y , t-Jq—ki 



(2) 

VqV = q • ((k - q)ncfc_ g + knc fe ). 



With q — k = |q — k| and k(i) = k — k • st as before. 
Cfe is the equilibrium direct correlation function connected to 
the structure factor via the Ornstein-Zernicke equation S k = 
1/(1 — nek) lf5lTl . From the expressions in dB16l ) one can now 
see the earlier proposed properties, Eq. ( lA14b . The function 
n q can schematically be written 

f°° 8 

n q (i) =aj 2 tf(t) + b dsj 2 s —g( 3 )h(t). (B17) 
Jo 8s 

The first term in (IB 17b corresponds to n q , the second term 

(2) (3) 

to n q and n q . f(t), g(t), h(t) are functions of tjr in the 
liquid and \ j\t in the glass. Eq. ( lA14b follows. The fact 
that the terms in ( IB17b start linearly with t and s respec- 
tively comes because in Eq. ( IB 16b is symmetric in k x , 
y{2) = v( 2 )(jfe2) ) and because at time t = (or s = 0) 
is anti-symmetric in k x , V^(—k x ) = —V^ l '(k x ), and the 
property n q (t — > 0) ~ 7 2 i follows after integration over 



<i 3 fc. The linear increase with time follows for example from 

^y(^) = 7^3;* 

For the numerical evaluation of Eq. (IB 16b , the transient 

correlator $ q (i) and the static structure factors S q and Sq 
are needed. As a purely technical simplification, we use the 
isotropic approximation ll36ll . which reads for long times in 
glassy states 

$ q (i) «*,(«) = f q e- c % m , (B18) 

with the non ergodicity parameter f q and the amplitude h q . 
The parameter c can be derived from a microscopic analysis, 
we use c = 3 IT36ll . For the static equilibrium structure factor, 
we use the Percus-Yevick closure 15111 . and approximate S q = 

Sq, which holds well at small shear rates fich . although the 
structure is nonanalytic ll80ll . In the limit of small shear rates, 
the contribution of the short time decay of the correlators to 
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the above expressions vanishes. The above expressions are 
evaluated using spherical coordinates with grid k max = 50, 
Afc = 0.05 and A9 = Acf) = tt/40 or smaller. The time grid 
in both t and s was jt = 2' 1 / 4 /10 10 , starting from i = 90 
corresponding to jt 6 10~ 4 . The results are included in 
Figs.E]and[T3] 



The right hand side of the vertex reads 



V 



s.2) 



qk 



qnc k . 



(B21) 



2. Incoherent Case 

From now on we denote incoherent functions with super- 
script s. The terms in Eq. (IA22b for the incoherent case, 
/ = Qa are approximated simi larly to the coherent analogs, 



using the pair density projector tl8 lh 



Pi(t) 



E 

P ,k 



£p(t)£k(t))(£p (t) £> k(t ) 



NS k 



(B19) 



k(t) 



The approximation for n q (t) is then straight forward fol- 
lowing Eqs. ( 1B2HB4I ). Regarding the vertex, there occur sim- 
plifications, 



V£ 1] = (v*yQ S *Q S Q*Q S ^) = k x k y ^ 



Si. 



(B20) 



For the memory functions n q (£) and 



>,3) 



the projector P| according to Eqs. dB7| ) and 



(£) we first use 
We note 



that the two terms corresponding to Eq. (IB lib vanish in this 
case independently. We arrive at the expressions equivalent to 
Eqs. dHE) and ( 15131 . reading (a xy exp[fit s ] e «*fitg«_ kek ) 
and (a xy exp[f2Ts](f2Tg»*)g»_ k k ^. We use the triple density 
projector as before, 



gk(t)gp(t)g°w)(gk(t)gp(t)g°W 
/ , ^77; 7 (pZZ) 

k>p>n 



TVS', 



«(*) 



according to Eq. (IB 15b . The discussions around Eqs. ( IB 16b 
and ( IB 16b and the approximations for the six-point functions 
hold similarly. We arrive at 
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I *k(t)*i_ q («), 

(B23b) 

(B23c) 



5", 



k • qn4. 



With = {g s q * g q ) / (nS q ) HH. We evaluate these expres- 
sions numerically, using the approximation ( IB 18b for the co- 
herent correlator and a similar approximation for the incoher- 
ent one, motivated by the solution for the correlator near the 
critical plateau J5(l HH]. We write for long times in glassy 
states 

Kit) « = /» cxp (-pc | 7 |t) , (B24) 

and again c = 3. We consider the case where the tagged 
particle has the same size as the bath particles, for which c 8 q = 
(S q - 1) /(nS q ) holds. We use the grid k rnax = 25, Afc = 0.2, 
AO = tt/40 and A</> = tt/20 (z-direction) and A6 = tt/20 
and A<fi = 7r/64 (y-direction) or smaller, the time is discretize 
as in the coherent case. 



The results are included in Figs.|2]and[T3]and again allow to 
conclude qualitative agreement with the approach of the main 
text, with quantitative differences arising from the involved 
MCT approximations. 

APPENDIX C: SIMULATION DETAILS 

The shear step model was simulated as follows: The one 
dimensional random walk within the well was discretized. At 
each time step (with At the length of the time step), the par- 
ticle position propagates by the step length s, x(t + At) = 
x{t) + s. s was determined by 

12 

s = J2 r i~6, (CI) 

t=l 
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with r,i random numbers (0,1). This gives a Gaussian distri- 
bution of width (7 = 1. We used a = 4 for the width of the 
well. If x(t + At) lies outside the well, the step is rejected, 
i.e., the particle stays at its position. For the case with exter- 



nal force F, the distribution in Eq. ( ICU was shifted by 5/100, 
corresponding to F = -^Z w j t h a 2 = 2D At. With a = 1 
and a = 4, the deviation of the linear response result is of the 
order of 1 percent. 



1] A. Einstein, Annalen der Physik 17, 459 (1905). 
2] R. Kubo, M. Toda, and N. Hashitsume, Statistitical Physics 2 

(Springer, Berlin, 1985). 
3] H. Nyquist, Phys. Rev. 32, 1 10 (1928). 
1] A. Crisanti and F. Ritort, J. Phys. A 36, R181 (2003). 
5] H. Risken, The Fokker-Planck Equation (Springer, Berlin, 
1984). 

5] G. S. Agarwal, Z. Physik 252, 25 (1972). 
7] T. Speck and U. Seifert, Europhys. Lett. 74, 391 (2006). 
S] V. Blickle, T. Speck, C. Lutz, U. Seifert, and C. Bechinger, 

Phys. Rev. Lett. 98, 210601 (2007). 
3] L. Berthier, J.-L. Barrat, and J. Kurchan, Phys. Rev. E 61, 5464 
(2000). 

3] L. Berthier and J.-L. Barrat, J. Chem. Phys. 116, 6228 (2002). 
1] L. Berthier and J.-L. Barrat, Phys. Rev. Lett. 89, 095702 (2002). 
2] J.-L. Barrat and L. Berthier, Phys. Rev. E 63, 012503 (2000). 
3] C. S. O'Hern, A. J. Liu, and S. R. Nagel, Phys. Rev. Lett. 93, 
165702 (2004). 

i] T. K. Haxton and A. J. Liu, Phys. Rev. Lett. 99, 195701 (2007). 
5] F Zamponi, G. Ruocco, and L. Angelani, Phys. Rev. E 71, 
020101 (2005). 

5] I. K. Ono, C. S. O'Hern, D. J. Durian, S. A. Langer, A. J. Liu, 

and S. R. Nagel, Phys. Rev. Lett. 89, 095703 (2002). 
7] J.-L. Barrat, J. Phys.: Condens. Matter 15, SI (2003). 
3] J.-L. Barrat and W. Kob, Europhys. Lett. 46, 637 (1999). 
3] W. Kob and J.-L. Barrat, Eur. Phys. J. B 13, 319 (1999). 
3] A. Latz, J. Phys.: Condens. Matter 12, 6353 (2000). 
1] P. Ilg and J.-L. Barrat, Europhys. Lett. 79, 26001 (2007). 
2] J. S. Langer and M. L. Manning, Phys. Rev. E 76, 056107 
(2007). 

3] N. Greinert, T. Wood, and P. Bartlett, Phys. Rev. Lett. 97, 
265702 (2006). 

i] B. Abou and F. Gallet, Phys. Rev. Lett. 93, 160603 (2004). 
5] C. Maggi, R. d. Leonardo, J. C. Dyre, and G. Ruocco, 

arXiv:0812.0740. 
5] C.GodrecheandJ. M.Luck, J. Phys. A: Math. Gen. 33, 1151 
(2000). 

7] C. Godreche and J. M. Luck, J. Phys. A: Math. Gen. 33, 9141 
(2000). 

3] P. Calabrese and A. Gambassi, Phys. Rev. E 65, 066120 (2002). 
?] P. Mayer, L. Berthier, J. P. Garrahan, and P. Sollich, Phys. Rev. 

E 68, 016116 (2003). 
3] P. Sollich, S. Fielding, and P. Mayer, J. Phys.: Condens. Matter 

14, 1683 (2002). 
1] F Corberi et al., J. Phys A: Math. Gen. 36, 4729 (2003). 
2] P. Calabrese and A. Gambassi, J. Stat. Mech.: Theo. Exp. p. 
P07013 (2004). 

3] P. Calabrese and A. Gambassi, J. Phys. A: Math. Gen. 38, R133 
(2005). 

4] M. Kriiger and M. Fuchs, Phys. Rev. Lett. 102, 135701 (2009). 
5] M. Fuchs and M. E. Cates, Phys. Rev. Lett. 89 (2002). 
5] M. Fuchs and M. E. Cates, Faraday Discuss. 123, 267 (2003). 
7] M. Fuchs and M. E. Cates, J. Phys.: Cond. Mat. 17, 1681 
(2005). 

[38] M. Fuchs and M. E. Cates, J. Rheol. 53, 957 (2009). 



[39] M. Fuchs, Advances in Polymer Science (2008), submitted, 

larXiv:08 10.2505 1 
[40] R. Besseling, E. R. Weeks, A. B. Schofield, and W. C. K. Poon, 

Phys. Rev. Lett. 99, 028301 (2007). 
[41] J. Zausch, J. Horbach, M. Laurati, S. Egelhaaf, J. M. Brader, 

Th. Voigtmann, and M. Fuchs, J. Phys.: Condens. Matter 20, 

404210 (2008). 

[42] F. Varnik. Complex Systems ed M. Tokuyama and I. Oppenheim 
(Amer. Inst, of Physics, 2008) p 160. 

[43] J. K. G. Dhont, An Introduction to Dynamics of Colloids (Else- 
vier science, Amsterdam, 1996). 

[44] I. Gazuz, A. M. Puertas, Th. Voigtmann, and M. Fuchs, Phys. 
Rev. Lett. 102, 248302 (2009). 

[45] P. Habdas, D. Schaar, A. C. Levitt, and E. R. Weeks, Europhys. 
Lett. 67, 477 (2004). 

[46] T. Harada and S. i. Sasa, Phys. Rev. Lett. 95, 130602 (2005). 

[47] E. D. Elrick, Austral. J. Phys. 15, 283 (1962). 

[48] G. Szamel, Phys. Rev. Lett. 93, 178301 (2004). 

[49] M. Kriiger, F. Weysser, J. Zausch, J. Horbach, and T. Voigt- 
mann, in preparation. 

[50] M. Kriiger, Properties of Non-Equilibrium States: 
Dense Colloidal Suspensions under Steady Shear- 
ing (PhD Thesis, Universitat Konstanz, 2009), URL 
[http : / / nbn- re solving . de/urn : nbn : de : bsz : 352-opus- 8 07 

[51] J. -P. Hansen and I. R. McDonald, Theory of Simple Liquids - 
2nd ed. (Academic press limited, London, 1986). 

[52] The violating term according to Eq. ((A3} reads Axq(i) = 
iV q (0)$ q (0)£ + Q (t 2 ). While $ q (0) = 1, it follows with 
Eqs. (TATJl and tBl6[ that 7V q (0) = OQj\). 

[53] F. Varnik and O. Henrich, Phys. Rev. B 73, 174209 (2006). 

[54] W. Gotze, Z. Phys. B 56, 139 (1984). 

[55] W. Gotze. Liquids, freezing and glass transition ed J.-P. Hansen, 
D. Levesque and J. Zinn-Justin (Amsterdam, 1991) p 287. 

[56] J. J. Crassous, M. Siebenbiirger, M. Ballauff, M. Drechsler, 
D. Hajnal, O. Henrich, and M. Fuchs, J. Chem. Phys. 128, 
204902 (2008). 

[57] M. Siebenbiirger, M. Fuchs, H. Winter, and M. Ballauff, J. 

Rheol. 53, 707 (2009). 
[58] More sophisticated versions of Eq. l[32} exist 03, [H], which 

are not necessary for our purpose. 
[59] F. Varnik, J. Chem. Phys. 125, 164514 (2006). 
[60] R. Graham, Z. Physik B - Cond. Mat. 40, 149 (1980). 
[61] J. A. McLennan, Introduction to Non-equilibrium Statistical 

Mechanics (Prentice Hall, New York, 1988). 
[62] M. Baiesi, C. Maes, and B. Wynants, Phys. Rev. Lett. 103, 

010602 (2009). 

[63] L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, 
Volume 6: Fluid Mechanics (Pergamon Press, Oxford, 1959). 

[64] D. J. Pine, J. P. Gollub, J. F. Brady, and A. M. Leshansky, Na- 
ture 438, 997 (2005). 

[65] J. P. Gollub and D. J. Pine, Physics Today 59, 8 (2006). 

[66] J. P. Bouchaud, J. Phys. 1 2, 1705 (1992). 

[67] A. Barrat and M. Mezard, J. Phys. I 5, 941 (1995). 

[68] C. Monthus and J. P. Bouchaud, J. Phys. A 29, 3847 (1996). 

[69] B. Rinn, P. Maass, and J. P. Bouchaud, Phys. Rev. Lett. 84, 5403 



21 



(2000). 

[70] S. Fielding and P. Sollich, Phys. Rev. Lett. 88, 050603 (2002). 
[71] S. M. Fielding, P. Sollich, and M. E. Cates, J. Rheol. 44, 323 
(2000). 

[72] C. W. Gardiner, Handbook of Stochastic Methods (Springer, 
Berlin, 1985). 

[73] We ignore the fact that the movement is discontinuous. 
[74] M. Kriiger and M. Fuchs, in preparation. 
[75] T. Speck and U. Seifert, Phys. Rev. E. 79, 040102(R) (2009). 
[76] W. Gotze and A. Latz, J. Phys.:Condens. Matter 1, 4169 (1989). 



[77] M. Fuchs and K. Kroy, J. Phys.:Condens. Matter 14, 9223 
(2002). 

[78] B. Cichocki and W. Hess, Physica A 141, 475 (1987). 
[79] M. Fuchs, W. Gotze, and M. R. Mayr, Phys. Rev. E 58, 3384 
(1998). 

[80] O. Henrich, O. Pfeifroth, and M. Fuchs, J. Phys: Condens. Mat- 
ter 19, 205132 (2007). 
[81] M. Kriiger et al, in preparation. 



